The condensation in non-growing complex networks under Boltzmann limit 
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We extend the Bianconi-Barabasi (B-B) fitness model to the non-growing complex network with 
fixed number of nodes and links. It is found that the statistical physics of this model makes it an 
appropriate representation of the Boltzmann statistics in the context of complex networks. The 
phase transition of this extended model is illustrated with numerical simulation and the correspond- 
ing "critical temperature" T c is identified. We note that the "non-condensation phase" in T > T c 
regime is different with "fit-get-rich" (FGR) phase of B-B model and that the connectivity degree 
distribution P(k) deviates from power-law distribution at given temperatures. 
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During the last decade there has been a burst of re- 
search activities on complex networks (see, e.g., [3, 0, 0, 
0, HI and references therein) , from the theoretical mod- 
eling as well as the empirical study of real-world net- 
works. Different approaches and models have already 
been proposed for growing complex networks @, 0, Hf, 
and the attention has been paid mainly to the out-of- 
cquilibrium dynamics by the means of preferential at- 



tachment @, M EL EE El- Of particular interest is 
the Bianconi-Barabasi (B-B) fitness model [l3| and the 
corresponding condensate behavior in the growing net- 
works (l3 |. 

In the fitness model, each nodes of the network corre- 
sponds to different energy levels of a system through the 
relation 



logr/j 



(1) 



where rji is the fitness of node i, and f3 is an inverse 
temperature parameter (i.e., /? = 1/T). The probability 
that a new link connects to the ith node depends on the 
fitness (energy level) rji and degree of connectivity (the 
number of links) of the ith node ki, 



n, 
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In Ref. 



141 ] . the evolving network with the fitness 
model was mapped into an equilibrium Bose gas. With 
mean-field arguments, the well-known Bose-Einstein 
(BE) statistics of the occupation number of an energy 
level (with energy e), n(e) = l/{e^ e ~ lx ' — 1), was ob- 
tained in the thermodynamic limit (t — > oo). As a re- 
sult, the network may undergo distinct phases - so called 
"scale-free" phase, "fit-get-rich" (FGR) phase and Bosc- 
Einstein-condensate (BEC) phase. The appearance of 
BEC phase indicates a "winner-takes-all" phenomenon - 
the fittest node has a finite fraction of the total number 
of links - during the evolution of the network. 

Although mapping to Bose gases successfully, one no- 
tices that the non-equilibrium (e.g., the growth of the 



number of nodes and the inertness of the links), and ir- 
reversible natures make B-B model an unreal representa- 
tion of BE statistics of Boson quantum gases. This mo- 
tivates us to extend the fitness model to a non-growing 
network with the fixed number of nodes and links. 

This non-growing network with fitness is constructed 
in the following way: initially, a graph composed of N 
nodes and K links (N > K) is given, and each node 
is assigned an energy chosen from some energy level 
distribution 17(e). We then let the network evolve: at 
each time step, a link is randomly picked up and discon- 
nected one end while the other end is kept fixed, then 
this link is rewired to a new node with the probability 
Hi (see eqn. J5])) until time steps reach the number of 
links K. The corresponding rate equation for the degree 
of connectivity of the ith node, ki, read as 
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where 
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is the partition function, and kio is the initial distribu- 
tion of connectivity degree at time t = 0. Eqn. is a 
group of ./V coupled differential equations, in principle, 
we can solve it and determine the degree of connectivity 
ki. However, a qualitative analysis is of interest here: 
first, the whole procedure can be seen as a mapping to 
Bose gases but with the characteristics of equilibrium 
evolution and the rewiring ability of links (so avoid the 
inertness in B-B model); second, in the low temperature 
regime, the rate equation can be simplified to 



dh _ fci(e ,t, fci )e-^ ei 
dt ~ Z t 



dkj_ 
dt 
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where k\ is the degree distribution of the lowest energy 
level. 

The form of eqn. ([5]) is similar to eqn. (3) in Ref. 14 1. 
However, it is noteworthy that there, eqn. (3) holds for all 
energy levels (and henceforth rusults in a BE statistics). 
Following eqn. ((SJ) and eqn. the links are eventually 
disconnected from higher energy levels, and rewired to 
the lowest energy level, this corresponds to a condensa- 
tion phenomenon 151 ] . and fci saturates when the two 
terms of right hand side (rhs.) of eqn. ((3]) cancels. 

The scenario of this non-growing network under the 
classical limit is very interesting and deserves detailed 
investigations. Under the classical limit, the links (or 
particles) arc independent of each other, and the connec- 
tivity probability of links depend only on the fitness rji or 
energy level e,; , but not on the degree of connection of the 
ith node, hi. We want to emphasize that this difference 
in the probability of connection makes the whole phase 
structure very different with FGR phase in the growing 
networks, as what we are going to explain later with our 
simulations. In the classical limit the connectivity prob- 
ability of network becomes 
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where 
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In classical limit, this non-growing network constructed 
in a given "temperature" (3 corresponds to an equilibrium 
gas composed of K independent particles distributed in 
N energy levels. According to the standard statistical 
physics jfor chemical potential [i and the partition func- 
tion Z [Ty] , we have 



-i In-. 

f3 K 



combine (jSJ and ([9|) we get 



Eqn. ([TO]) simply leads to 

* = £*< = £< 



(9) 



(10) 
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which is, clearly, Boltzmann statistics with distribution 

h~ f{ei) = e-^-»\ (12) 

In order to show the phase transition of this non- 
growing network, we numerically simulate the above sce- 
nario under the classical limit, we generate a graph con- 
sist of N nodes and K links. These links arc rewired 
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FIG. 1: Simulation results of the network under Boltzmann 
limit for distribution 17(e) = Ce B with — 0. The upper 
panel is the plot of absolute value of chemical potential fi 
with temperature T on a logarithmic scale, the lower panel 
is the fraction of total links occupied by the most connected 
nodes. The curves are plotted for the same number of nodes 
N — 1000, but different number of total links: K = 500 (open 
triangle with solid line), K = 200 (open square with dashed 
line), and K — 100 (open circle with dotted line), respectively. 
The numerical simulations are averaged over 50 runs. 



to existing nodes by the following way: one point of ev- 
ery link connect to one node j randomly while the other 
point connect to another node i with probability IT by 
eqn. ([7]). The situation that two or more links connecting 
the same pair of nodes is not allowed. 

The simulation results of our model in the classical 
limit are shown in Figs. Q]- El In Fig.[IJ we show the re- 
lationship of absolute value of chemical potential \fi\ with 
temperature parameter T = 1//3, in upper panel. The 
curves represent networks with N = 1000 and K = 500, 
200, 100, respectively. It can be seen clearly that there 
exist transitions in which means the chemical poten- 
tial /x changes its sign, i.e., /1 is positive in T <T C regime 
and negative in T > T c regime. The negativeness of /z in 
T > T c regime only has a formal similarity to that of the 
ideal Boltzmann gases, since, after all, the classical limit 
condition is not satisfied in the low temperature regime 
and hence no such condensation behavior in Boltzmann 
statistics. 
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FIG. 2: Similar to the Fig.[T] but for the distribution of energy 
level g(e) = Ce e with 9 = 1. 



The (non)condensation behavior is shown in the lower 
panel in Fig. [TJ We can see that the fraction increases 
to unity gradually, as temperature decreases. Also, the 
critical temperature T c shifts towards to higher temper- 
ature when the ratio K/N increases. In low temperature 
limit, the connectivity probability of non-growing net- 
work is primarily dominated by factor e~^ £i rather than 
ki, which is similar to a network under the classical limit. 
As we already mentioned earlier, it should be noted that 
the non- condensation phase of classical limit network is 
different from the so called "FGR" phase of the fitness 
model in growing network. In non- condensation phase 
the connectivity probability of each node is randomly dis- 
tributed over all energy levels. However, in FGR phase, 
for the nodes with higher fitness the connectivity degree 
increase more quickly than those nodes with lower fitness. 

Different distributions of energy level, g{e), also in- 
fluence the (non)condensation phase behavior under the 
classical limit. In our model, we employ the following 
form, 

5 (e) = Ce 9 , (13) 

where 9 is a parameter and C is normalization factor. 
The simulations shown in Fig. [T] correspond to 9 = 0, 
C = 1, i.e., a uniform distribution. In Fig. [5J we show the 
simulation results of and k max /K for = 1, C = 2. 



FIG. 3: The comparisons of the absolute value of chemical 
potential, (upper panel), and k max /K (lower panel) with 
the same number of nodes (N = 1000) and links (K = 200), 
but with different energy level distributions g(e) = Ce° , 9 = 
(open diamond with solid lines), and 9 = 1 (star with dashed 
lines). 



To illustrate clearly the effects of g(e), in Fig.[3]we show 
the comparisons of the simulations with 9 = and 9 = 1 
for the same number of nodes and links (N = 1000, K = 
200). The enhancement of critical transition temperature 
T c for distribution with 9 = 1 is obvious. 

We also plot the probability distribution of connectiv- 
ity degree P(k) at temperature T = 0.2 for distributions 
with 9 = (Fig.Q} and 9=1 (Fig.©, respectively. The 
up left, up right and lower left panels in Fig. Q] and Fig.O 
correspond to the following number of nodes and links: 
N = 1000 and K = 500, 200, 100. With different num- 
ber of nodes N and links K but the same ratios of K/N 
as shown in the up left and lower right panels in Fig. [4] 
and Fig. El respectively, the distributions P{k) (vs. k) 
are quite similar except for the fluctuation in tails. A 
particularly interesting observation is, at given temper- 
ature, it seems like the probability distributions of fc; 
deviate from power-law distribution, as with the increase 
of the ratios K/N , which is clearly demonstrated in both 
figures. This significant difference from scale free net- 
work show some hints on the possible connection with 
recent observations: some real- world complex networks 
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FIG. 4: The connectivity degree distribution P(k) vs. degree 
k for N — 1000, K = 500 (open triangles), 200 (open squares), 
100 (open circles) and N = 2000, K = 1000 (open diamonds) 
with distribution g(e) = Ce e , 9 = at temperature T — 0.2. 





FIG. 5: Similar to the Fig.|4]but for the distribution of energy 
levels g(e) = Ce e with 9 = 1. 



display a Weibull- or mixed- Weibull-power-law distribu- 
tion 0, EE EE 0, ED, [H]. Although some models are 



proposed to explain (see, e.g., Ref. [l8|), the underly- 
ing mechanism for such distribution is worthy of further 
studying. 

In conclusion, in present paper, we extend Bianconi- 
Barabasi fitness model from growing complex networks 
to non-growing networks with fixed number of nodes and 



links. We fulfill numerical simulations and identify the 
corresponding "critical temperature" T c , We also find 
that the "non-condensation phase" of our extension dif- 
fers from FGR phase in B-B model, and the degree dis- 
tribution P{k) deviates from the power-law distribution 
at some temperatures. 
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